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' The study of the relative phase between strong and electromagnetic amplitudes is of great 

importance for understanding the dynamics of charmonium decays. The information of 
the phase can be obtained model-independently by fitting the scan data of some special 
decay channels, one of which is tp' ~ > K+ K~ . To find out the optimal data taking 
strategy for a scan experiment in the measurement of the phase in if)' — > K + K~, the 
minimization process is analyzed from a theoretical point of view. The result indicates 
that for one parameter fit, only one data taking point in the vicinity of a resonance peak 
is sufficient to acquire the optimal precision. Numerical results are obtained by fitting 
simulated scan data. Besides the results related to the relative phase between strong and 
electromagnetic amplitudes, the method is extended to analyze the fits of other resonant 
parameters, such as the mass and the total decay width of ip' . 

Keywords: e + e — annihilation, relative phase, statistical optimization 

PACS numbers: 02.60.Pn, 13.25.Gv, 13.40.Hq 

1. Introduction 

The charmonium hadronic decay is mainly through two processes: the strong and 
the electromagnetic interactions. The relative phase between the strong and the 
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electromagnetic decay amplitudes is an important parameter in understanding de- 
cay dynamics. Studies ha ve be en carried out for many J/tp tw o-body decay modes: 
Vector-Pseudoscalar (VP) 0^1, Pseudoscalar-Pseudoscalar (pp)EHEI Vector- Vector 
(VV) El and Nucleon-antiNucleon (AW) El. These analyses reveal that there exists 
a relative ortho gonal phase between the strong and the electromagnetic amplitudes 
in J/tp decays l l | 2 | 3 | 4 | 5 | 6 | 7 [ ^ ^ ^here is also a th eoretical argument which 
favors the ±90° phase El. Experimentally, some analyses 1 9 1 1 | 1 1 | based on limited 
VP and PP data indicate that such a phase is compatible with the data. More- 
over, some efforts have b een made to extend the phase study to V(3770) decay 
phenomenologically l l^ | 13 | anc j experimentally ^31. 

The phase study can provide valuable clue for exploring the relation between 
the stro ng a nd the electromagnetic interactions. Now with the upgraded accelerator 
BEPCIlQSl and detector BESIE^ a luminosity of 6.5 x 10 32 cm" 2 s" 1 has achieved, 
which is the highest luminosity in r-charm energy region ever existed. 22 6 M J/tp 
events, 106 M ip' events, and 2.9 fb^ 1 ^(3770) data have been collected even 
more colossal data are to be collected in the forthcoming years, which gives a great 
opportunity to determine the phase between the strong and the electromagnetic 
amplitudes with unprecedented statistical precision. 

However, examining the existing determination of the relative phase, since the 
data are merely taken at one or two energy points, we find most of studies are 
model-dependent. A typical model assumption is the SU(3) symmetry in charmo- 
nium decays which supply additional constraint on the electromagnetic decay am- 
plitudes in charmonium decays into similar final states such as VP, PP and so on. 
Now with a high luminosity accelerator, it is possible to measure the phase model- 
independently by scanning the cross sections in the vicinity of the resonance. As 
the strength of the resonance decays varies with energy, the precision of the phase 
measurement depends on the data taking energy when the total data taking time 
is fixed. Therefore, the optimization study for the data taking strategy is of great 
importance in order to obtain the most precise results with the limited luminosity 
(equivalently within the limited data taking time). 

Without losing generality, we focus on the mode of ip' decays to K + K~ final 
state. Because, as will be shown in the next section, this decay mode can accom- 
modate a comparatively simple parametrization form which is of great benefit to 
extract the relative phase. 

As far as the optimization of data taking strategy is concerned, sampling sim- 
ulation technique was adopted for optimizing the r mass measurement I 18 | 19 [ ^n 
interesting conclusion from the study is that for one parameter fit, data at only one 
energy point is enough to acquire the best precision. The breakthrough of this mono- 
graph lies in that the minimization process is analyzed in detail from a theoretical 
point of view, which leads to the same conclusion as that of r mass measurement. 
As a cross check, numerical results are obtained by fitting simulated ip 1 scan data. 
Moreover, this method is extended to extract other resonance parameters, such as 
the mass and the total decay width of ip' . 
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2. Theoretical Framework 

For xjj' — > 0~0~, the tt + tt~ channel is through electromagnetic decays, the KgKl 

through SU(3) breaking strong decays, and the K + K~ through both. Therefore, the 

— > K + K~ decay is the only process which can be used to study the phase between 

strong and electromagnetic interactions in an energy scan experiment. Taken into 

ac count th e continuum process, the decay amplitude of this mode is parametrized 
191101201 

A K+K - =E C + E+^M 1 (1) 

where E c is the continuum amplitude, E the electromagnetic amplitude, and n^-M 
the SU(3) breaking strong amplitude. They can be expressed explicitly as 

E c ot-, Eoc-B(s), ^-M <xCe**-B(s), (2) 
s s 2 s 

where the real parameters <f> and C are the relative phase and the relative strength 

between the strong and the electromagnetic amplitudes, and B(s) is defined as ^ 

= ynyq ( } 

w s - M 2 + iMT t ' v ' 

where %fs is the center of mass energy, a is the QED fine structure constant, M 
and T t are the mass and total width of ip', T ee is the partial width of ip' — » e + e~. 
The Born order cross section for this channel reads 

°k° + "k- = [1 + mC^B{s)) + |C B( S )| 2 ] 

x\F K +k-(s)\ 2 Pk+k-{s), 
where = 1 + Ce'*; F K+K -{s) is the form factor, which is usually written as 

_2_ 

'■is 



f k+k-( s ) = Jk+k-I s with Ik+k- being a constant; P k +k-( s ) = it Ik 1S tlic 
phase space factor, with 



2 jp2 2 " 2 

Ik = e k -m K =-- ™ K , 

where qx is the momentum of K + or K~ , uik is the nominal mass of K meson. 

In actual experiment, the effect of Initial State Radiation (ISR) is considered 
through an integral 

r x f 

oy.c(s) = / dxF(x,s)<7 B om(s(l - x)), (5) 



where Fix, s) is the structure function which can be calculated to an accuracy of 

o,i% M 
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In addition, another important experimental effect, the energy spread of e + and 
e~ must also be taken into co nsideration. Finally, the experimentally observed cross 
section is expressed as | 24 | 25 | 

/•OO 

Vex P (Vs) = / dv / s 7 cr r . c .(v / s 7 )G(\/s 7 , v 7 ^), (6) 
Jo 

where G(v?, y/s) is a Gaussian distribution 

r- r- 1 (v^-y^) 2 

G(vV, y/s) = -^^e 2A^. 
V2^fA 



Here A indicates the energy spread of the collision beams. 

Some parameter values for the numerical calculation in the following sections 
are articulated in Table [T] 

Table 1. Some parameter values for 
numerical calculation. The quantity 
with * will be set as a free fitting pa- 
rameter for the corresponding study. 



Quantity numerical value Remark 

*M 3.68609 GeV Ref. ED 

*r< 304 keV Ref. EH 

r ee 2.35 keV Ref. ED 

m K 493.677 MeV Ref. ED 

f K+K - 0.9 GeV 2 Ref. 23 

A 1.3 MeV Ref. EH 

±<f> 90° Ref. 123 

C 2.5 Ref. 123 



3. Minimization Analysis 

For a scan experiment, several points, say totally N p t points, need to be taken 
in a vicinity o f a resonance (in this monograph the ip'). The estimator is usually 
constructed as ^5 



N, 

,2 _ 



Y 2 =y 1 1 l ' (7) 

' (AN obs ) 2 ' 



where N° bs and AN° bs = yj N° bs arc the observed number of events and its error 
at the i-th point, Li the corresponding luminosity, £i the selection efficiency, and 0$ 
the theoretical cross section that is a exp in Eq. ©. The fitting parameters (relative 
phase, strength, etc.) are contained in <7j, and these parameters and the correspond- 
ing errors can be extracted by minimizing the x 2 function defined in Eq. ([?])• In the 
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following analyses, only concerned is one free fitting parameter, the relative phase 
between strong and electromagnetic amplitudes, that is, cj>. 

If we denote the observed cross section measured at energy point i as a° bs , and 
rewrite 



N. 



obs 



t -obs 



Eq. (JTJ) can be recast as 



/ = E SS« bs - "J 2 = ^EJ« is - « 



of 



(8) 



(9) 



Here % 2 is replaced with / for simplicity, and the following relations are utilized : 



XiL , 



1 



(10) 



where Lq is the total luminosity (corresponding to the finite total data taking time) 
and Xi is the fraction of luminosity at the i-th energy point. Moreover, Si is supposed 
to be the same at all points (e = 50% is used for numerical calculation), which is a 
fairly good approximation for the scan of narrow resonances, such as J/ip and tp' . 

In the light of Eq. (j9|), the first and second order derivatives of the function / 
to <b can be derived as 



df 



obs 



(11) 



dcj) 2 



2L Q e 



a obs 



doi 



I „obs \ 



d 2 a t 



(12) 



The experimentally concerned cross section functions are generally smooth 
enough, which can be approximated by polynomial functions. Therefore, the first 
and second order derivatives of these functions are also smooth enough (refer to 
Fig. [1] when the parameters take the values in Table [TJ. Under such case, we argue 
that the second term in Eq. (fl~2|) could be neglected. When the fitting process fin- 
ishes, <7i in Eq. (fTTj) and Eq. (fl2"|) can be considered as the true value of the cross 
section at energy point i. As we assumed previously, a° bs is the experimentally 
measured cross section at energy point i, then {a° bs — a{) could be considered as a 
random variable which satisfy a Gaussian distribution with mean as and deviation 
as Aof bs (the error of cr° bs ). As a conservative estimation, we assume the relative 
error of cross section measurement is 10%, which means Aa° bs = 0.1 • o° bs . The 
expectation of the second term inside the sum in Eq. ([T2"|) could be calculated by 
using sampling method. The comparison of these two terms are shown in Fig. [TJ 
from which the second term is quite small compared with the first one, therefore its 
effect can be neglected. Now Eq. (fT2|) becomes 
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(a) 



(b) 



Fig. 1. (a) Comparison between the cross section and its first and second order derivative. The 
solid line shows the cross section (<r). The dashed line shows the first order derivative (cr'). The 
dottcd-dashed line shows the second order derivative (cr"). (b) Comparison between first and 
second term inside the bracket in Eq. 1121 . The solid line only shows the value of the first term 
(Ti); the dotted line contains contributions both from the first term and the second term (Ti + 
T2). Values in Table [l] are taken in making these plots. 



21 



d(jj 
d<j> 



(13) 



where </>* is the fitting result of the relative phase and a* is the theoretical cross 
section at <p = 4>* ■ 

The fitting error of <f> can be evaluated as 



(14) 



According to Eq. (|T4|) , the maximum of the second order derivative of fitting func- 
tion yields the minimum of fitting error. Define a new function g as 

2 



1 {den 

9i = < lw 



(15) 



where the subscript i denotes the value of g at the z-th energy point. Then Eq. ([T 
becomes 



21 



= 2L e^ 



(16) 



Notice that Yli ac* = 1, it is readily to obtain the following inequalities 

grain ^ ^i)Qm.vn — ^ , %i9i — ( ^ •Ej)9max — 9maxi (i'O 

i i i 

where g m in (gmax) is the minimum (maximum) value of g within the energy region 
concerned. To get maximum J^. a^fft, only one data taking point is sufficient and it 
should be located at the energy point which renders g maximum. 
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4. Numerical results 

To reinforce the preceding conclusion, the simulated scan data are fit to get numer- 
ical results. In this procedure, great many times of fitting need to be performed, 
where the large number of calculations must be carried out for the observed cross 
section. Unfortunately, two nested integrations of the observed cross section, which 
take into account both the ISR correction and beam energy spread effect, take 
so much time that any actual optimization fitting becomes impractical. In a re- 
cent study SH, using the generalized linear regression approach, a complex energy- 
dependent factor is approximated by a linear function of energy. Taken advantage 
of this simplification, the integration of ISR correction can be performed and an 
analytical expression with accuracy at the level of 1% is obtained. Then, the orig- 
inal two-fold integral is simplified into a one-fold integral, which reduces the total 
computing time by two orders of magnitude. In the following studies, the simplified 
observed cross section formulas are adopted to acquire all numerical results. 



4.1. Relative Phase 

Considering the parameter to be analyzed is the relative phase between strong and 
electromagnetic amplitude of tf>' decay, the distribution of g^ and the fitting error 
Ecf, on energy region when 4>* = 90° is shown in Fig. [21 according to which, at the 
energy point 3.6868 GeV, the value of function g^ reaches its maximum while 
reaches its minimum. In the vicinity of 3.684 GeV, the g^ value is very small and 
the corresponding E$ is quite large. So this point (3.684 GeV) should be avoided 
in the scan experiment 

By fixing the energy point to 3.6868 GeV, the error obtained from fitting and 
computed by Eq. (|T4"|) is shown in Fig. [3] Just as expected, the higher the luminosity, 
the smaller the error. Moreover, the fitting and computing values of error are so 
consist with each other that it is hardly to distinguish them in Fig. [3] (a) . To exhibit 
the details, the relative difference of E^,, that is, 

_ ^(computing) - ^(fitting) 
^(computing) 

is shown in Fig. [3] (b). 

If the relative phase <f> variates, the optimal position of energy will change cor- 
respondingly. Table [5] lists the optimal values of energy position for some special 
phase angles. According to these information, the values of the optimal energy points 
arrange from 3.686 GeV ~ 3.687 GeV, nearly within the scope of 1 McV. 



a To validate this result, the sampling technique is used to check the data taking point distribution 
and the fitting error. Details about the sampling technique can be found in Appendix. 
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Fig. 2. The distributions of gj, and for <j>* = 90°. 




Fig. 3. (a) Comparison of the equivalent error for different luminosity by fixing the energy point 
at 3.6868 GeV; (b) Relative difference R<p at 3.6868 GeV for different luminosity. 



Table 2. The optimal data 
taking position for different 
relative phase </> 

<t>(°) Optimal point (GeV) 



3.68604 

45 3.68700 

90 3.68680 

135 3.68706 

180 3.68648 

270 3.68672 
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4.2. Other parameters 

The minimization analysis discussed in Section [3] is applicable to any parameter 
we are concerned with. As long as the variable <f> is replaced with the parameter 
to be analyzed, all aforementioned deductions are valid. In study that followed, we 
perform the optimization for the two interested resonance parameters, mass and 
total width. 



4.2.1. ip' Mass 

The similar analyses are performed for the mass of ijj' resonance, and results are 
displayed in Fig. [6j where four curves corresponding to different relative phases, 
4> = 0°, 90°, 180°, and 270°. There are two new features for the optimization of 
mass parameter. Firstly, the energy position for smallest Em is at 3.6845 GeV for <p 
= 0°, 90°, and 270°; but at 3.6874 GeV for = 180°. Secondly, two energy positions 
should be avoided due to larger values of Em- One is around 3.686 GeV, the other 
is near 3.68 GeV for = 0°, 90°, and 270°; while near 3.69 GeV for = 180°. 

To understand the heterogeneous behavior of curve for <f> — 180° from the other 
ones, we take the curve for <fi = 90° as a representative, and show the first order 
derivatives of cross sections in Fig. [5j where the left and right rows correspond to 
90° and 180°, respectively. The total cross section a exp is divided into three parts, 
i.e. c^p, cr^p, and cr^ p , as we did in Ref. ^1 The corresponding derivatives of 
them are shown sequently in Fig. [5] Investigation of those figures indicates that the 
crucial role for the different behavior between 90°- and 180°-curves is played by the 
crf^p, the derivative variation of which is opposite to each other. Furthermore, if we 
scrutinize the equations in Section 3 of Ref. the sign of the derivative of a^ p 
is determined by the coefficient = 6(r ee /a) • (1 + C cos0), which is positive 
when 4> = 90° and negative when <fi = 180° for C — 2.5. The switching point is at 
(f> = cos _1 (-l/C) w 113.5°. 
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A=l.lMcV - 

A=1.3MeV - 
A=1.5MeV - 



A=l.lMeV - 

A=1.3MeV - 
A=1.5MeV - 



(a) 



(b) 



Fig. 6. The results of analyzing tp' total width for A = 1.1 MeV, 1.3 MeV and 1.5 MeV. (a) is 
the gr value and (b) is the fitting error Er- 



4.2.2. ijj 1 Total Width 

The optimization results for the total width of ip' resonance are shown in 
Fig. [51 where three curves corresponding to distinctive beam energy spreads, 
A = 1.1,1.3,1.5 MeV, are presented. With the enhancement of A, the position 
of minimum error, Ep, shifts a little bit rightward along the abscissa. For all cir- 
cumstances, the energy position for the maximum Er is near 3.68 GeV while for 
the minimum error is around 3.686 GeV. 



5. Conclusion and Discussion 

In this paper, one-parameter-and-one-point conclusion is demonstrated through a 
theoretical analysis of minimization process instead of the sampling simulation as 
we did before. As far as the phase study is concerned, for the ifj' —> K + K~ process, 
the optimal data taking point is determined to at 3.6868 GeV which is near the 
peak of ip' nominal mass. The same method is also used to acquire the optimal 
point for other resonance parameters, such as the mass and the total width of ip' . 

In principle, the idea put forth in Section [3] could be extended for multi- 
parameter optimization. Formally speaking, the vector and matrix quantities, would 
be adopted for the corresponding analysis, say, the second order derivative of one 
parameter are to be replaced by Hessian matrix, a matrix of second order derivative 
for all parameters. 

However, there are some problems not easily to be settled. The most prominent 
one is how to define "optimal" . In one parameter scenario, the optimal data taking 
point is the one which could make the fitting error of the parameter reaches its 
minimum. But for multi-parameters, there are many options: the sum of relative 
fitting errors of all parameters reaches its minimum; the merely fitting error of one 
major parameter reaches its minimum while others do not. Different options lead 
to distinctive results. All this makes the situation more complicated and is left to 
the study in the further. 
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Appendix A. Sampling Technique Methodology 

Suppose there are N pt data taking points in experiment, and the theoretical number 
of events in i-th energy point could be calculated as 

Nt e = Li ■ en ■ e,i = 1,2,..., N pU (A.l) 

where Li is the integrated luminosity in i-th energy point, Ei is the event selection 
efficiency 

In sampling technique, the experimentally observed number of events and its 
error could be taken as 

AN° bs = (N? 16 ) 1 / 2 , (A.2) 

N° hs = Nl he + £ ■ AN° bs , (A.3) 

where £ is a random number which satisfy Gaussian distribution. 

Using the observed event number and its error calculated above, the parameter 
we interested in (in this paper, the relative phase) could be fitted by finding the 
minimum of Eq. (J7J . 

By repeating this process, a large number of observed event number and error 
could be generated and so does the fitting parameter and its error. We can compare 
these errors with the results obtained from the method we just developed. 

In section we scan through the energy region using one data taking point. 
To check this result, five consecutive energy points are used in energy scan. The 
fitting error of five consecutive points scan versus the central energy point is shown 
in Fig. [7] This result is similar with one point scan. Their difference is also shown in 
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Fig. [7J The difference between the two scan schemes is generally at the level of one 
per mille except for the points around 3.684 GeV, where the variation of error curve 
is rather rapidly. As to the five-point scheme, there is at least one point within the 
region with comparatively large error. 



References 

1. DMII Collab. (J. Jousset et al), Phys. Rev. D 41, 1389 (1990). 

2. Mark III Collab. (D. Coffman et al), Phys. Rev. D 38, 2695 (1988). 

3. M. Suzuki, Phys. Rev. D 60, 051501 (1999). 

4. G. Lopez, M.J.L. Lucio and J. Pestieau, |hep-ph/99 02300 

5. L. Kopke and N. Wermes, Phys. Rep. 174, 67 (1989). 

6. R. Baldini et al, Phys. Lett. B 444, 111 (1998). 

7. M. Suzuki, Phys. Rev. D 63, 054021 (2001). 

8. J.M. Gerard and J. Weyers, Phys. Lett. B 462, 324 (1999). 

9. C.Z. Yuan, P. Wang and X.H. Mo, Phys. Lett. B 567, 73 (2003). 

10. P. Wang, C.Z. Yuan and X.H. Mo, Phys. Rev. D 69, 057502 (2004). 

11. BES Collab. (J.Z. Bai et al), Phys. Rev. Lett. 91, 052001 (2004). 

12. P. Wang, C.Z. Yuan and X.H. Mo, Phys. Lett. B 574, 41 (2003). 

13. P. Wang, X.H. Mo and C.Z. Yuan, Int. J. Mod. Phys. A 21, 5163 (2006). 

14. BESII Collab. (M. Ablikim et al), Phys. Rev. D 70, 077101 (2004). 

15. Q.Qin et al, Status and Performance of BEPCII, IPAC-2010-WEXMH01, May 2010. 

16. Lhe BESIII Collaboration, Design and Construction of the BESIII Detector, arXiv: 
0911.4960, 2009. 

17. ZY. Wang, G. Li, K. Zhu, et al, ip' total number, BAM-00003; H.X. Yang, B.X. 
Zhang, X.T. Liao, et al, Determination of J/ip total number, BAM-00011; Lili Jiang, et 
al, Measurements of the luminosity at 3.773 GeV and 3.650 GeV, BAM-00037. 

18. Y.K. Wang, X.H. Mo, C.Z. Yuan, et al, Nucl. Instr. Meth. A 583, 479 (2007). 

19. Y.K. Wang, JY. Zhang, X.H. Mo, C.Z. Yuan, et al, Chinese Physics C 33, 501 (2009). 

20. P. Wang, X.H. Mo and C.Z. Yuan, Phys. Lett. B 557, 192 (2003). 

21. E.A. Kuraev and V.S. Fadin, Sov. J. Nucl. Phys. 41, 466-472 (1985). 

22. G. Altarelli and G. Martinelli, CERN 86-02, 47 (1986). 

23. F.A. Berends, G. Burgers and W.L. Neerven, Nucl. Phys. B 297, 429 (1988); F.A. 
Berends, G. Burgers and W.L. Neerven, Nucl. Phys. B 304, 921 (1988). 

24. SY. Lee, Accelerator Phystcs, 2nd edn. (FuDan University Press, Shanghai, 2006) 

25. K. Wille, The Physics of Particle Accelerators (Oxford University Press, New York, 
2000) 

26. Particle Data Group (K. Nakamura et al), J. Phys. G 37, 075021 (2010). 

27. CLEO Collab. (T.K. Pedlar et al), Phys. Rev. Lett. 95, 261803 (2005). 

28. BES Collab. (J.Z. Bai et al), Phys. Lett. B 550, 24-32 (2002). 

29. CLEO Collab. (S.Dobbs et al), Phys. Rev. D 74, 011105 (2006). 

30. Y.S. Zhu, Probability and Statistics in Experimental Phystcs, 2nd edn. (Science Press, 
Beijing, 2006) 

31. B.Q. Wang, X.H. Mo, et al, Chinese Physics C 35, 411 (2011). 



